MZ-TH/96-30 
hep-ph/9610285 
September 1996 



New representation of the two-loop crossed vertex function 



Alexander Frink, Ulrich Kilian, Dirk KreimeiQ 

Institut fur Physik, Johannes Gutenberg- Universitat, 
D-55099 Mainz, Germany 



Abstract 

We calculate the two-loop vertex function for the crossed topology, and for 
arbitrary masses and external momenta. We derive a double integral represen- 
tation, suitable for a numerical evaluation by a Gaussian quadrature. Real and 
imaginary parts of the diagram can be calculated separately. 
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1 Introduction 



In the recent past many efforts have been made to calculate two-loop Feynman dia- 
grams with masses, e.g. [jl], |2|, IS, ji, EH. Though now very effective methods exist for 
the general mass case of the two-loop two point function || [7J, these methods need 
considerable extensions or modifications to cope with two-loop three-point functions 
in general. The planar topology has been discussed extensively, e.g. in [|, |3|, ]8|, ||. 
For the other important topology — the crossed topology — so far two methods have 
been presented. 

The method presented in |J is based on Feynman parameters and uses high 
dimensional Monte Carlo integration, resulting in extensive CPU usage and slow 
convergence when high accuracy is needed. 

Taylor expansion, analytic continuation and Pade approximations, as presented 
in 0, gives results with very high accuracy, but the method is, at this stage, still 
restricted to one kinematical variable. 

Our approach is based on @, [§, where it gave excellent results for the planar 
two-loop three point function. The calculation of the crossed topology is similar 
to the planar case. Again we succeed by considering orthogonal and parallel space 
variables. From this starting point, a four- fold integration is immediately obtained. 
Still following the lines of || we use Euler transformations for subsequent integra- 
tions. The difference between the planar and the non-planar topology is apparent in 
the increasing number of different cases which have to be considered for the latter, 
while, fortunately, the conceptual frame remains unaffected. Once more we end up 
with a two-fold integral over a finite region, solely involving dilogarithms and related 
functions. This integral representation allows for a Gaussian quadrature, and is thus 
perfectly suitable for practical purposes |H| . 



2 Calculation 




Figure 1: The two-loop crossed vertex function. 
Our task is to calculate the two-loop crossed vertex function, as shown in fig.|l]. 
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We restrict ourselves to the scalar integral, which is given by 

where q = q\ + q 2 is the momentum of the incoming particle and q\ and q2 are the 
momenta of the outgoing particles. Further, I and k are the internal loop momenta, 
and Pi, i = 1 . . .6 are the propagators of the internal particles, as labeled in fig.[l]. 
We assume that all external particles are massive and have time-like momenta. The 
massless limit is approached smoothly as long as no on-shell singularities occur. 

The integral will be splitted into parallel and orthogonal space integrations. The 
parallel space is two-dimensional here, hence leaving two dimensions for the orthog- 
onal space. The integral is convergent in D = 4 dimensions, and there is no need to 
use a regularization scheme. 

Since the integral is Lorentz invariant, the splitting into parallel and orthogonal 
space variables is equivalent to the choice of a Lorentz frame. For time-like external 
particles, we choose the rest frame of the incoming particle, and assume that outgoing 
particles are moving along the x-axis. 

Explicitly, the four-momenta are chosen as 

<? = (ei + e 2 , 0,0,0) 
Qi = (ei,fc,0,0) 

<?2 = (e 2 ,-g 2 ,o,o) 

l» = {lo,hJ±) 

h" = (k ,k u k ± ) (2) 

with q z , e\ and e 2 positive and q z < e\, q z < e%. 

The propagators Pi of the internal particles can be written down using the flow 
of loop momenta as shown in fig.|I] 

Pi = (I + k — qi) 2 — m\ + it] 

P 2 = (l + k + q 2 f -mj + iri 

Ps = {I - qi) 2 - m\ + 17] 

Pa = (k + q 2 ) 2 -m\ + irj 

P5 = I 2 — m 2 + ir] 

P 6 = k 2 -m 2 6 + iri. (3) 

A small imaginary part — irj with rj > is assigned to the squared masses of the inter- 
nal particles and is chosen equal for all particles for convenience. Hence differences of 
propagators have a vanishing imaginary part, which we utilize in the following. The 
choice of different (small, positive) 77, 's for the propagator would not change our final 
result, but is technically more cumbersome. 

Using the explicit representation of the momenta in eq.(^), a shift of the loop 
momenta with a trivial Jacobian according to 

h -> h + q z 
h — > Iq + h 

k -> ko + h, (4) 
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achieves that the propagators are linear in l\ and k\ and hence can be written as 



Pi = (Z + k - df + 2(/ + k - ei) (Zj + fci) -l\-k\- 2l ± k ± z -m 2 x + ir) 

P 2 = {l + h + e 2 f + 2(l + h + e 2 ) (h + h) - l\ - k\ - 2l ± k ± z - m\ + ir] 

Ps = {l - erf + 2(l - ei )k - l 2 ± - m 2 3 + ir] 

P 4 = {k + e 2 + q z f + 2(k + e 2 + q z )(h - q z ) - k]_ - m\ + vq 

P5 = (l ~ q z ) 2 + 2(l - q z )(h + q z ) -ll-ml + tr] 

P 6 = fcg + 2/cqA;! - k\ - ml + i-q , (5) 

where z is the cosine of the angle between l± and k±. 

The volume element in the integral in eq.(P is given by (cf. 0) 

1 dz 
d I d A k = —dlo dko dl\ dk\ ds dt da j= (6) 



z 



with s = l±, t = k\ and the angle a describing the absolute position of l± and k±. 
The integration over a gives a trivial factor 2-7T. 

An important difference between the crossed and the planar two-loop vertex func- 
tion is that now necessarily two (instead of one) propagators depend on z, which, in 
our notation, are Pi and P 2 . This is a result of the fact that in the planar case the 
loop momenta can be arranged to flow through only one common propagator, which 
is not possible for the crossed topology. 

But after applying a partial fraction decomposition to the integrand 

11 1 

+ 77^ F7T77777777^ ( 7 ) 



P 1 P 2 P 3 P i P 5 P 6 (P 2 - Pi)PiP 3 P 4 P5P 6 (Pi - P 2 )P 2 P 3 P 4 P 5 P 6 

the ^-dependence in P x — P 2 drops out and the ^-integration can be performed as in 
the planar case: 

1 

dz 1 7T 

k = l,2 (8) 



i VI - z 2 A k + B k z + ir] ojA 2 k -B 2 k + 2A k i v 



where Pi = A\ + B\z + ir] for the first term in eq. (0) and P 2 = A 2 + B 2 z + it] in 
the second, respectively. <fz denotes the square root of a complex number with a cut 
along the positive real axis, whereas yfx is the usual square root of a positive real 
number (if x has a small imaginary part of any sign, it can be ignored). 

The next step is the integration over Zi and k\ using Cauchy's theorem for both of 
them. It is important to notice that, as a result of the partial fraction decomposition 
in eq. ([?p, in the difference Pi — P 2 the imaginary part of the masses —irj, which has 
been chosen to be equal for all masses, drops out. Therefore not all poles in l\ and 
k± lie in the upper or lower complex half-plane, but some also on the real axis. The 
integral can be made meaningful if interpreted as a principal value integral .0 Then 

2 One can avoid principal value integrals by choosing different imaginary parts r\i for the propa- 
gators, as mentioned earlier. After having convinced ourselves that the results remain unchanged, 
we prefer to follow the route outlined here, for purely technical reasons. 
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one has to use the modified Cauchy's theorem [[11]]: 



oo 



n 



P.V. / f(z) dz 



2 ?7 r£Res(/(z)) 



z=Zi,lm(zi) sC 



— oo 



+ivr]TRes(/(z)) 



z=Zj, lm(zj)=0 



(9) 



j 



where the z% are the n poles of f(z) inside the closed integration contour, whereas Zj 
are m poles along the integration path. Poles on the path contribute only with half 
the weight compared to poles inside the path. For the cases we are confronted with 
it is guaranteed that the function tends to zero sufficiently fast for large l\ and ki, 
hence the theorem is applicable. 

The integration path has to be closed either in the upper or in the lower complex 
half-plane, depending on the position of the cuts of the square root $ A 2 — B 2 + 2Air\. 
Similarly to the planar case, the position of the cuts is determined by the sign of 
+ ko — ei for the first term in eq.(0) and by the sign of lo + k + e 2 for the second 
term. 

Let us concentrate on the first term of the partial fraction decomposition, the 
second can be handled analogously. Assume l + k — 1\ < 0, so that we have to close 
the contour of the l\ integration in the lower half-plane. Then, some propagators 
have poles inside the path: P 3 iff l — e\ > and P5 iff lo — q z > 0. Additionally, 
P 2 — Pi always contributes with half its residue, since it has a pole on the real 
axis. In the following we will call residues, which involve only P3 . . . P@, complex 
contributions, because all poles lie in the upper or lower complex half-plane, whereas 
residues involving Pi — P 2 are called real contributions, because the pole lies on the 
real axis. 

For the ki integration, we close the contour in the same half-plane as for the 
li integration, depending on the sign of Iq + k — e%. This is not necessary for the 
P 2 — Pi contribution, where the square root becomes independent of ki, but it is done 
for consistency. Let us first assume that we had a Zi-pole in P 3 . Then we have poles 
in ki from P4 iff ko + e 2 + q z > and P§ iff ko > 0, and a pole in P 2 — Pi , iff lo — &i < 0. 
But the last case is a contradiction to lo — t\ > above, so it does not contribute. 

Analogously, if we had a /i-pole in P5, we have the same constraints from P4 and 
Pq, but lo — q z < from P 2 — Pi, which is again a contradiction. 

The last /i-pole in P 2 — Pi, gives contributions from all other propagators, 

namely P 3 iff l — ei < 0, P 4 iff k + e 2 + q z > 0, P 5 iff l - q z < and P 6 iff k > 0. 

If l + k Q — ei > 0, the relation operators have to be reversed in all inequalities. Out 
of the complex combinations, only the pairs (P 3 ,P 4 ), (P 5 ,P 4 ) and (P 5 ,P 6 ) contribute 
in triangular regions in the (ko, lo) plane if lo + k — e\ < (see fig.0). If lo + k — e% > 0, 
no terms contribute. All real residues contribute for lo + ko — ei < as well as for 
lo + ko — ei > 0. The areas are not triangles, but unbound, as can be seen in fig.^. 

The second term of the partial fraction decomposition also gives us three con- 
tributing complex residues (all for lo + k + e 2 > 0, fig.(4]) and four real residues, which 
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Figure 2: Triangles from the first term of the partial fraction decomposition [complex 
contributions) . 



have the same constraints as for the first term, but with l Q + k — e% ^ replaced by 

k + k + e 2 § (fig.||). 

If we compare the real contributions from both terms of the partial fraction de- 
composition eq.(M), we notice that the terms for two corresponding residues (e.g. 
(P 2 — Pi,P$) an d (-Pi — P21P3)) are equal, except for an overall minus sign. As a 
consequence they cancel in the area where both contribute together. This is every- 
where the case except in the strip — e 2 < Iq + k < e\ (see fig.§). Furthermore it 
can be shown that all four real contributions from fig|| add up to zero where they 
all contribute together. This is the case everywhere outside the finite area shown in 
fig.|7|, which is also the joined area of all triangles in figsfj and |] from the complex 
residues. 

After the residue integrations we end up with a four-fold integral representation 
which is similar in nature, but slightly more complicated in its technical appearance 
than the planar case. It can be written in the general form 



00 00 

V(gi,92) = Yl JJ di odk jds fdtC(k ,l 



3 Aj 

1 1 



+ 61 + cis a 2 t + b 2 + c 2 s a 3 t + 63 + c 3 s ^( a t + b + cs) 2 - Ast 



(10) 



The j-sum runs over all areas Aj depicted in figs.|2| - §. The coefficients C(k , l ), di, 
bi, Ci (i — 1 . . . 3), a, b and c are in general dependent on j, l and k . bi and b have 
an infinitesimal imaginary part, all other coefficients are real. Some coefficients are 
vanishing: for the complex residues, d\ and £3 are zero, and we have one pure pole 
in s, one in t and a mixed (s,t) pole. For the real residues, either c 3 = or a 3 = 0, 
i.e. we have two mixed poles and either one pure pole in t or in s. In the latter case 
we can relabel s<->£, together with exchanging di <-» q and a <->■ c, so that we always 
have exactly two poles in s, either pure or mixed, and C3 = 0. In the planar case we 
always had one pure s pole and two pure t poles, with mixed poles altogether absent. 
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Figure 3: Unbound areas from the first term of the partial fraction decomposition 
(real contributions). 

To proceed further with the s and t integrations, we apply a partial fraction 
decomposition in s, followed by a similar partial fraction decomposition in t. 

1 1 1 



with 



ait + bi + CiS a 2 t + b 2 + c 2 s a 3 t + b 



a [-2 1 ) 1 !__\ av 

J + t 02 t + t J \s + s 02 {t) S + Soi{t)J 



c 

SQl(t) 



1 



c 2 (ai&3 - 61P3) - ci(a 2 6 3 - b 2 c 3 ) 
h 

C 2 6i - Ci& 2 



c 2 ai — cia 2 
ait + 6i 

Cl 




(P> . PJ (Pa ■ Pe) (P, ■ Pe) 



Figure 4: Triangles from the second term of the partial fraction decomposition (com- 
plex contributions). 

_ a 2 t + b 2 , . 

S 2^) = — = • (12) 

c 2 

Thus we have to calculate four integrals of the form 

00,00 

v=r^L[ * _ 1 da) 

J q t + t J s + s {t) ^ ( a t + b + cs) 2 - 4st 

with So(t) being either a linear function of t or a constant. Now we have to split these 
into real and imaginary part. It can be shown that b always has a positive small 
imaginary part, whereas the sign of the imaginary part of so and to is a function of 
lo and ko, so that in general sq = so ± it] and to = to ±in. In contrast, in the planar 
case the imaginary part was always negative. 

Now there are three possible sources for an imaginary part of V: either s > ^0 
or the argument of the square root can become negative. In the first two cases the 
imaginary part can be extracted using 



LXJ (JXJ (JXJ 

lim n / T ~T ■ fa) = P-V. / T ™ I S(x + x )f(x)dx . (14) 
v->o J x + xo ± vn J x + xq J 
' 



To analyze the contribution from the square root, we have to distinguish between 
complex and real residues. For the complex case, the area in the s — t plane, where 
the argument is negative, is an ellipse. For a given positive t, the s values on the 
border of the ellipse can be calculated as 



1 

a V2 = ~2 



Vi±Jt(l - ac) - be 



(15) 



which has real and positive solutions for < t < bc/(l — ac) if b > 0, since 1 — ac < 0, 
cf. fig. H For the real residues, a = — 1 and c = — 1, therefore the area where the 
argument becomes negative is unbounded with a parabola as its boundary (fig. |9|). 



S 



1 1 





Figure 5: Unbound areas from the second term of the partial fraction decomposition 
(real contributions). 



Now we can do the s-integration with the aid of an Euler's change of variables 
[ |H| and obtain 



V 



t + tf 



(/(s , <r») + ik(t)g(s , 0*)) dt 



(16) 



with 



and 



f(s ,ai 



7^ arctan ( £2 ygp) if 0i < -s < a 2 
i_ soifcyo-i^+ v^ 



else 



if (Ti < -s < cr 2 



A else 



00 



W 
(ii) 



(17) 



(18) 



where we note that So and o~i are functions oft and further R = Sq + SQ(a\ + a 2 ) + 0-^2 
and k = —2, —1,0, 1 or 2, depending on the relative magnitudes of so and a\/ 2 and 
the sign of the imaginary part of s . One has to choose + A /o"iO'2 in eqs.([L7D and ( |l"8"D 
if t > —b/a, and —^fd^a^ else. The condition o\ < — s < cr 2 is equivalent to R < 0. 
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Figure 6: Sum of both terms of the pfd for the real contributions (dark shading from 
the first, light shading from the second term = — first term.) 



The ^-intervals where k is constant, together with the corresponding k values, can 
be calculated by solving so = and so = — <?i for t. With another Euler's change of 
variables, using eq.(P^) and exploiting the relevant properties of the functions log(x) 
and arctan(a;), the t-integration leads to expressions of the form 



X2 



P.V. / -f- — from g(s , a t ) , (19) 

x 2 + px + q 



XI 

X2 



f arctanirx + s) „. . ... . . 

P.V. / - 1 ; from / (s Q ,<Ti ), case i 20 
J x 2 + vx + a 



and 



PV. f lQ gl x +rx + s\ 



x 2 + px + q 



which can be expressed in terms of logarithms, arcus-tangens, dilogarithms and 
Clausen's functions [19, as it was the case for the planar topology. The full re- 



sult is, as expected, a lengthy expression in terms of these special functions, which 
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Figure 7: Sum of all four terms from fig. |6|. 



b>0 




Figure 8: The ellipse. 

we cannot list here. We rather follow the philosophy of || and present examples in 
the next section. 

We stress that to obtain stable numerical results in the two dimensional integral 
over l Q and k it is very important to add all contributions (real and complex) to the 
same Z and k , because eq.(^) introduces some artificial divergences which cancel in 
the sum. 

As expected, the numerical integrand for the crossed topology is of the same 
nature as for the planar case, but involves more terms and different cases. This 
naturally increases the amount of CPU time needed to obtain the requested accuracy: 
as a thumb rule, we found that the crossed topology demands 5-10 times more time 
compared to the planar case. 

The threshold behaviour can be examined with the four-fold integral representa- 
tion eq . (|ToD . The crossed vertex function has three two-particle thresholds and six 
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Figure 9: The parabola. 



three-particle thresholds. As stated above a possible source for an imaginary part is 
a negative argument of the square root. A necessary condition for this to happen is 
b > for some Iq, ko inside the integration region. It can be shown that each three- 
particle threshold corresponds to one of the six b coefficients of the complex triangles 
|J. The other b coefficients belonging to the real regions are identical and correspond 
to the two-particle threshold q 2 > (mi+mj) 2 . The remaining two-particle thresholds 
q\ > (m 3 + m 5 ) 2 and q\ > (m 4 + m 6 ) 2 can be identified with the coefficients s i 
(when di = 0) and t i respectively which produce an imaginary part when s i and 
toi become negative. A possible four-particle threshold q 2 > (m 3 + m 4 + m 5 + mo) 2 
disconnects the graph into three parts and is therefore a combination of the two 
two-particle thresholds q\ > (m 3 + m 5 ) 2 and q\ > (m 4 + me) 2 . 



3 Examples 

Only few analytical and/or numerical results are known for the crossed vertex function 
to compare with. Besides the symmetries with respect to internal masses and external 
momenta, some limiting cases can be checked. In the case of zero momentum transfer 
the crossed vertex function reduces to the master two-point topology |], [7| with a 



squared propagator. fig.|10] shows these limits for some arbitrary masses (mi = 2, 
m 2 = 3, m 3 = 1, m 4 = 1.5, = 0.6 and m 6 = 0.2), q 2 = 1 and either qi — >• or 
q2 —>■ 0. Additionally, for vanishing q 2 we obtain the vacuum bubble calculated in 

For the case of all internal masses vanishing an analytical formula in known |14| . 



This limit is approached smoothly in fig.[ll] for all masses going to zero simultaneously 
and ei = 4, e 2 = 3 and q z = 1. 

In the limit of the momenta of the outgoing particles on the light-cone (qf = 
q\ = 0) and all internal masses equal this diagram can be calculated with the small 
momentum expansion technique [||. A comparison for real and imaginary part is 
show in fig.0. We can obtain this limit easily by transforming our parallel space 
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Figure 10: Ratio of crossed vertex- function to master two-point function for qi (short 
dash) and q 2 — > (long dash). 
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Figure 11: Limit of all masses vanishing. 

coordinates to light cone coordinates, which does not interfere with our subsequent 
steps. 

As a last example in fig.[TB| we show the decay Z — ► tt with the exchange of two 
Z bosons. This diagram has been calculated in by a five- dimensional numerical 
integration over the Feynman parameters. 



Conclusions 

In this paper, we demonstrated the calculation of the two-loop vertex function for the 
crossed topology. We outlined how to achieve the reduction to a manageable two-fold 
integral representation, and verified its correctness by comparison with the literature, 



wherever data were available. Further, our integral representation was used in flO 
with results in perfect agreement with the expectations. 
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Figure 12: All masses equal, q\ = q\ = 0, real part (left) and imaginary part (right). 
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Figure 13: Decay Z — > tt, real part. 



Our results complement the results for the planar topology in ||. Together, 
the scalar two-loop three-point function is now available in D = 4 dimensions for all 
topologies, all masses, and arbitrary external momenta. Some degenerated topologies 



were already given in |]10| , obtained by similar methods. Such degenerated cases 
typically also appear when one confronts tensor integrals. For the future, we plan to 
incorporate such cases in the package XLOOPS |J. Also, code which implements the 
results presented here will be incorporated there. 

To our knowledge, no other method is at this stage able to deliver reliable results 
for the massive two-loop vertex function in such generality and accuracy. 
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